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Transverse momentum spectra up to 4.5GeV/c around midrapidity of tt*, p, p, in Pb+Pb 
reactions were measured at ^s„„ — 17.3 GeV by the CERN-NA49 experiment. The nuclear mod- 
ification factors Raa for tt^ and Rcp for ,p,p, were extracted and are compared to RHIC 
results at ^/T^ = 200 GeV. The modification factor Raa shows a rapid increase with transverse 
momentum in the covered region. This indicates that the Cronin effect is the dominating effect in 
our energy range. The modification factor Rcp, in which the contribution of the Cronin effect is 
reduced, shows a saturation well below unity in the tt* channel. The extracted Rcp values follow 
the 200 GeV RHIC results closely in the available transverse momentum range, except for vr* above 
2.5 GeV/c transverse momentum. There the measured suppression is smaller than that observed at 
RHIC. 
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I. INTRODUCTION 



One of the most interesting features discovered at 
RHIC is the suppression of particle production at high 
transverse momenta in central nucleus-nucleus reactions 
relative to peripheral ones as well as to p+nucleus and to 
p+p collisions |E S Hi S Hi ■ This is generally interpreted 
as a sign of parton energy loss in hot and dense strongly 
interacting matter. 

The aim of the presented analysis is to investigate the 
energy dependence of these effects via a systematic study 
of Pb+Pb reactions at top ion-SPS energy, 158v4GeV 
(y^;;^ = 17.3 GeV), with the CERN-NA49 detector. A 
similar study has been published by the CERN-WA98 
collaboration for the 7r° channel The aim of this 
paper is to extend their results to all charged particle 
channels, i.e. tt^, p, p and . 

Invariant yields were extracted as a function of trans- 
verse momentum in the range from 0.3 to 4.5GeV/c 
in the rapidity interval —0.3 <y< 0.7 (midrapidity) , at 
different collision centralities. The identification of par- 
ticle types is crucial, because the particle composition of 
hadron spectra changes rapidly with transverse momen- 
tum and differs significantly from that observed at RHIC 
energies (at SPS energies p production is comparable to 
the TT^ yields already at moderate values oi p^ ). 

The most important results are the nuclear modifica- 
tion factors RaAi which were extracted from the iden- 
tified charged hadron spectra in Pb-fPb and from the 
already published p-|-p pion spectra Q- It is important 
to note that the use of reference spectra (e.g. p-|-p), mea- 
sured exactly at the same collision energy is necessary, 
because of the rapid change of the shape of the parti- 
cle spectra around y/s^^ = 17.3 GeV. A baseline, which 
is constructed from the existing nearby energy measure- 
ments is not sufficient, as discussed in Q. 

Besides the in-medium parton energy loss, the Raa 
quantity can also contain other medium effects, such as 
the Cronin effect a term used for the observation of 
an increased particle yield at high transverse momentum 
in p-|-nucleus compared to the elementary p-l-p reactions. 
This effect is believed to be a consequence of multiple 
scattering of projectile partons (or hadrons) in the nuclei 
before the particle production process. The Cronin effect 
could cause a large part of the nuclear modification effects 
measured by the Raa quantity. For a more sensitive 
search of the other nuclear effects the Rcp (central to 
peripheral) modification factors were also extracted from 
the data. If the multiple scattering interpretation of the 
Cronin effect is correct, its contribution to the Rcp ratio 
is expected to be reduced, as it should already be present 
in the peripheral baseline. 



II. EXPERIMENTAL SETUP 

The NA49 detector [l^ is a wide-acceptance hadron 
spectrometer for the study of hadron production in col- 
lisions of hadrons or heavy ions at the CERN-SPS. The 
main components are four large-volume Time Projection 
Chambers, TPCs, (Fig. [T]) which are capable of detect- 
ing 80% of some 1500 charged particles created in a cen- 
tral Pb-|-Pb collision at 158A GeV beam energy. Two 
chambers, the Vertex TPCs (VTPC-1 and VTPC-2), 
are located in the magnetic field of two superconduct- 
ing dipole magnets (1.5 and 1.1 T), while the two others 
(MTPC-L and MTPC-R) arc positioned downstream of 
the magnets symmetrically to the beam line. The setup 
is supplemented by two Time of Flight (TOF) detector 
arrays, which are not used in this analysis, and a set 
of calorimeters. The NA49 TPCs allow precise mea- 
surements of particle momenta p with a resolution of 
a(p)/p2 (0.3 - 7) • 10-4 (GeV/c)-i. 

Pb foils (typically of 224mg/cm^ thickness) are used 
as target for Pb-|-Pb collisions, and a liquid hydrogen 
cylinder (length 20 cm) for p-l-p interactions. The target 
is positioned about 80 cm upstream from VTPC-1. 

Pb beam particles are identified by means of their 
charge as seen by a Helium Gas-Cerenkov counter (S2') 
and proton beam particles by a 2 mm thick scintillator 

(52) . Both detectors are situated in front of the tar- 
get. For p beams, interactions in the target are selected 
by requiring a valid incoming beam particle and no sig- 
nal from a small scintillation counter (S4) placed on the 
beam line between the two vertex magnets. For p-l-p 
interactions at 158 GeV this counter selects a (trigger) 
cross section of 28.5 mb out of 31.6 mb of the total in- 
elastic cross section. For Pb-ion beams an interaction 
trigger is provided by a Helium Gas-Cerenkov counter 

(53) directly behind the target. The S3 counter is used 
to select minimum-bias collisions by requiring a reduc- 
tion of the Cerenkov signal. Since the Cerenkov signal 
is proportional to Z^, this requirement ensures that the 
Pb projectile has interacted with a minimal constraint on 
the type of interaction. This setup limits the triggers on 
non-target interactions to rare beam-gas collisions, the 
fraction of which proved to be small after cuts, even in 
the case of peripheral Pb-|-Pb collisions. The resulting 
minimum-bias trigger cross section was about 80% of the 
total inelastic cross section crinci = 7.15 b. 

For Pb-|-Pb reactions, the centrality of a collision is se- 
lected (on-line for central Pb-|-Pb, off-line for minimum- 
bias Pb+Pb interactions) by a trigger using information 
from a downstream calorimeter (VCAL), which measures 
the energy of the projectile spectator nucleons [ill ■ 



III. ANALYSIS 

TableUhsts the statistics used in the 158A GeV Pb+Pb 
collision analysis. 
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FIG. 1: (Color online) Setup of the CERN-NA49 experiment showing different beam definitions and target arrangements. 



Reaction 


[GeV] 


Centrality 


Number of events 


Pb+Pb 


17.3 


central 


830 k 


Pb+Pb 


17.3 


mid-central 


1.4 M 


Pb+Pb 


17.3 


peripheral 


200 k 



TABLE L The event statistics, used in the analysis. 



A. Event Selection 

As the target setup of the experiment is not contained 
in a vacuum pipe, non-target background collisions are 
also expected to be recorded in the unfiltered event sam- 
ple. To reduce this contamination, a cut was applied 
to the reconstructed longitudinal position of the colli- 
sion point: only those collisions were accepted, which 
occurred close to the nominal target position. An exam- 
ple plot is shown in Fig. [2] for the case of Pb-|-Pb reac- 
tions, showing the target peak and the small flat non- 
target contribution, which was used to estimate the re- 
maining background events. The contamination of the 
peripheral Pb-t-Pb spectra remains below 5%. For non- 
peripheral spectra this contamination is zero. The pe- 
ripheral Pb-fPb spectra are corrected for this effect. 

In the case of Pb-|-Pb reactions, the recorded events 
were classified by collision centrality, which is correlated 
to the impact parameter - the transverse distance of the 
centers of the colliding nuclei at impact. The centrality 
is determined via the energy measured by the VCAL: the 
non-interacting (spectator) part of the projectile nucleus 
travels along the beamline, finally hitting the VCAL's 
surface, and leaving an energy signal in the apparatus, 
proportional to the volume of the projectile spectator 
matter. The higher this energy is, the more peripheral is 
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FIG. 2: (Color online) The distribution of the longitudinal co- 
ordinate of the reconstructed collision point in minimum-bias 
Pb-f Pb events, together with the cuts. The other peaks cor- 
respond to material in the beamline, while the small constant 
background under the target peak corresponds to beam-|-gas 
reactions. 



the given reaction. The event centrality is defined by the 
fraction of total inelastic cross section, i.e. by the running 
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integral 

a(£vcAL)_ 1 /■^"'^"^ ME'yCAl) 

^ / , T^' Cl-t^VCAL; 

C^Incl CTlncl Jo "i-^VCAL 

where ainei is the total inelastic cross section, and — 
is the differential cross section in E^cAh , which is defined 
by the i?vcAL spectrum, normalizing its area to the trig- 
ger cross section. The VCAL energy spectrum and event 
centrality classification is shown in Fig. [3l A relatively 
large centrality interval was used for the selection of pe- 
ripheral events in order to gain enough statistics in the 
particle spectrum analysis. 
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FIG. 3: (Color online) The distribution of the energy de- 
posited in the VCAL by the projectile spectators in Pb-fPb 
minimum-bias reactions. The centrality intervals selected by 
fraction of total inelastic cross section are shown by the ver- 
tical lines. 

Due to the aging of the photomultiplier tubes of the 
calorimeter, the VCAL response was slowly varying in 
time in terms of scale offset and amplification. This ef- 
fect is corrected by a time-dependence calibration pro- 
cedure, which is based on keeping the correlation of the 
VCAL energy to the total measured charged multiplicity 
time-independent (see [l^). The mean values of vari- 
ous collision parameters (number of wounded nucleons, 
number of binary collisions) for the selected centrality 
intervals were calculated from the VENUS-4.12 Monte 
Carlo model [l3| which was filtered by the detector sim- 
ulation. For the non-peripheral region, this was straight- 
forward by applying statistical averaging in the centrality 
intervals. For the peripheral region, a semi-empiric ap- 
proach was used, which also takes the trigger bias into 
account: the VCAL energy distribution for the averaging 
was taken from the measurement, while the mapping of 



VCAL energy into collision parameters was taken from 
VENUS -I- detector simulation. The procedure is de- 
scribed in [1^. A list of the average values is shown 
in Table [IT] together with their systematic errors, which 
were estimated by assuming a 5% uncertainty of the to- 
tal inelastic cross section, and a 10% uncertainty of the 
skin thickness of the nuclear density profile. The stated 
values were also confirmed by the GLISSANDO Monte 
Carlo model |l3| ■ 



Bin 


Centrality 


{Nw) 




1. 


0-5% 


357 ±2 


742 ± 18 


2. 


5 - 12.5% 


288 ±3 


565 ± 18 


3. 


12.5 - 23.5% 


211 ±4 


379 ± 16 


4. 


23.5 - 33.5% 


146 ±5 


234 ± 14 


5. 


33.5 - 43.5% 


87±7* 


120 ± 13 * 


6. 


43.5 - 80% 


40±4 * 


46±5 * 


5.U6. 


33.5 - 80% 


56±7* 


70 ± 13 * 



TABLE II: The mean values of the number of wounded nu- 
cleons and binary collisions in various centrality intervals in 
Pb+Pb collisions, together with their systematic errors. *: 
Averages for the full minimum-bias dataset using the semi- 
empiric method, as discussed in the text. 



B. Track Selection 

The transverse momentum spectra of particles in 
heavy-ion collisions are known to decrease rapidly (ap- 
proximately exponentially) toward higher transverse mo- 
menta. Therefore, the background of fake tracks becomes 
more and more important with increasing , especially 
for fixed-target experiments, in which the track density is 
strongly focused in the forward direction, increasing the 
probability for fake track formation. We apply cuts to 
minimize the contribution of such artificial tracks. These 
background tracks consist of discontinuous track candi- 
dates, and track candidates whose trajectories pass close 
to the border of the detector volume. The former kind 
of tracks mainly arise from erroneous pairing of straight 
track pieces in MTPC, outside of the magnetic field, to 
some residual points in the VTPCs. The latter kind of 
tracks are subject to distortions, as they are likely to 
have missing or displaced clusters. The corresponding 
trajectories will have wrong fitted curvature, therefore 
they should be discarded from the analysis. 

The discontinuous track candidates can be easily iden- 
tified: they are track candidates that (according to their 
reconstructed charge and momentum) would have left 
hits in a given TPC (VTPCl, VTPC2, or MTPC), but 
did not. These tracks were rejected as one source of fake 
tracks. After this cleaning procedure a 3 dimensional 
phase space study was performed. Its coordinates are ra- 
pidity y, transverse momentum and charge-reflected 
azimuth (j). Charge-reflected azimuth is deflned by the 
azimuth angle of the momentum vector in the transverse 
plane, whose x-component is multiplied by the charge 



sign of the track candidate - the a;-direction being per- 
pendicular to the magnetic field and to the beamline. 
As our detector has an a:-reflection symmetry, the use of 
charge-reflected azimuth enables us to distinguish tracks 
not crossing the plane defined by the beamline and the 
magnetic field direction (right-side tracks) from tracks 
crossing this plane (wrong-side tracks). An example of 
the track quality distribution as a function of (j) and is 
shown in the upper panel of Fig. [4] for a typical rapidity 
slice. A domain of high quality right-side tracks reveals 
itself: in this region, the fraction of those tracks with 
a ratio of the number of measured to number of poten- 
tial points below 60% is very low. Selection of this re- 
gion by a 3 dimensional momentum space cut provides a 
clean track sample. The 3 dimensional momentum space 
cut is guided by the isosurfaces of the number of poten- 
tial points, and by the requirement of avoiding efficiency 
holes, which can be determined from the dropping of the 
(j) distribution at fixed {v^Pt) values (the 4> distribution 
should be uniform due to the axial symmetry of particle 
production). The momentum space cut curve is shown 
by the dotted line in the upper panel of Fig. IH at given 
y. The distribution is shown in the lower panel of Fig. 
ID at given {y,p^), which reveals the efficiency holes to 
be avoided, and the high track yield excess at the accep- 
tance borders. The dotted arrows indicate the momen- 
tum space cut in 0, for this particular (y^ p ) s lice. The 
outlined procedure is discussed in detail in ]l7| . The re- 
sulting phase space, which contains the clean track sam- 
ple, covers the rapidity region —0.3 < y < 0.7, the trans- 
verse momentum region OGeV/c < Pj. < 5GeV/c, and 
a {y,p^) dependent <j) interval. In the following analysis, 
all distributions will be shown as a function of transverse 
momentum. 

The particle identification was performed via the spe- 
cific energy loss (^) of the particles. In each trans- 
verse momentum bin the ^ spectrum was recorded, and 
by using the known shape of the response function (see 

[3 Ha]) 

of the given particle species {-k^ , p,p, K^, e^), 
a fit was performed to the ^ histogram, with the ampli- 
tudes of the particle response functions as free parame- 
ters. Also other parameters characteristic of the response 
function shape, such as the most probable value, were 
kept as free parameters. Thus, their statistical errors 
are also reflected in the fltted errors of the amplitudes. 
As our data was limited by statistics at higher p^, , Pois- 
son maximum-likelihood fitting was employed. This ap- 
proach allows reliable fits even in regions with low statis- 
tics (see \2d\). An example can be seen in Fig. [5] which 
demonstrates the ^ fit procedure for inclusive particle 
identification in a typical phase space bin. The full pro- 
cedure is described in [l7|. Systematic errors, caused 
by possible systematic shifts of fit parameters were es- 
timated by calculating the statistical regression matrix 
from the fitted covariance matrix, and by propagating the 
systematic errors of the fit parameters (approximately 
estimated to be about 0.5%) to the fitted yields via this 
matrix. The resulting systematic errors of yields are es- 
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FIG. 4: (Color online) Upper panel: example plot for the mo- 
mentum space selection in the rapidity slice < j/ < 0.1, with 
TT* mass hypothesis. (f> means charge-reflected azimuth. The 
size of the boxes indicate the fraction of those track candi- 
dates in the given momentum space bin, which have a ratio 
of number of measured to number of potential points below 
60% (i.e. the fraction of those track candidates, which most 
likely do not correspond to real particle trajectories). The 
color map indicates the number of potential points as a func- 
tion of momentum space. As can be seen, the bad track can- 
didates mainly populate the borders of the acceptance. The 
dotted line shows the momentum space cut, guided by poten- 
tial point isosurfaces. Lower panel: example plot for the <j> 
distribution at < y < 0.1 and 2.5 GeV/c < < 2.6 GeV/c. 



timated to be 4% for K+ , 2% for p, and 1% for 7r+, tt", 
p, K'. 
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FIG. 5: (Color online) Example plot for the identification of 
negative particles by specific energy loss fits in central Pb+Pb 
collisions (0-5%). The parameterization of the response func- 
tion of each particle species with fixed momentum is known, 
but the yields (the amplitudes in the summed ^ spectrum 
model) are fit parameters. 



C. Corrections 

Monte Carlo methods including a GEANT simulation 
of our detector were used to correct the yields of 7r+, tt^ 
and K~^,K~ as well as p and p for tracking inefficien- 
cies (below 10%), decay losses (from 20% to 0%), feed- 
down from weakly decaying particles (5 — 30%), and fi- 
nally for geometric acceptance. The peripheral particle 
spectra were also corrected for non-target contamination 
(5%) by extracting the non-target contribution from pure 
beam-|-gas events, and subtracting it (for detailed proce- 
dure, see [l3|). The yield of fake tracks, after the in- 
troduced track cuts, turned out to be negligible. The 
momentum space resolution was estimated to be better 
than 1% in the whole region. The momentum scale un- 
certainty is of the order of 0.1%. 

The largest and most sensitive correction is the feed- 
down correction. The yield of feed-down particles was 
estimated by using a simulation to determine the condi- 
tional probabilities of reconstructing a secondary particle 
originating from a weakly decaying primary particle as a 
primary particle. To calculate the feed-down yields as a 
function of momentum, one has to fold these tables of 
conditional probabilities with the momentum differenti- 
ated yields of the weakly decaying particles. The rele- 
vant decay channels are listed in Table IIIII The yields of 
the weakly decaying particles K^, A, A were taken from 
parameterizations of yields, measured previously in our 



experiment [Ml, [H, ^ {K° yield was constructed from 
average charged Kaon data) . The used A, A yields in- 
clude the S° — ^ A7 and the E*^ A7 contributions, but 
are feed-down corrected for S (S) decays. Since the S 
(S) yields are generally below 20% of the A (A) yields, 
their contribution is neglected in the correction. How- 
ever, it constitutes a source of systematic errors, that 
is discussed below. The K'^ yield does not suffer from 
feed-down contamination. 
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< 30% to p 
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p 7r+ 


< 30% to p 


E+ 


^ Ptt" 


< 6% to p 


E+ 


n TT^ 


negligible 


E" 


— * n TT~ 


negligible 


E- 




< 6% to p 


E" 


— > n TT~ 


negligible 


E+ 


n -R^ 


negligible 



TABLE III: The list of relevant feed-down channels. 

For 7r+, only the tt+tt^ channel gives a sizable 

contribution, while for n~ also the A — > pTr~ channel has 
to be taken into account. All other contributions to the 
TT* channels are negligible. For p, the A pir^ gives the 
dominant contribution, and for p, the A — > pi:^ is domi- 
nant. The p and p spectra are also contaminated by the 
decays of and E~. These are taken into account by 
scaling up the A, A yields by 20%. This treatment was 
suggested by the VENUS-4.12 model, which predicts that 
the relative intensity of the non-A contribution to the p 
feed-down, or the relative intensity of the non-A contri- 
bution to the p feed-down is approximately constant at 
20%. Assuming a 50% systematic uncertainty of this 20% 
scaling, the contribution to the systematic uncertainty of 
the p,p spectra would be 3%. There is an additional un- 
certainty of the p feed-down correction, which is due to 
the large errors of the measured A inverse slope parame- 
ters. The systematic error, caused by the poor knowledge 
of the A slopes was estimated by repeating the feed-down 
calculation using the A temperatures, which should pro- 
vide a reasonable upper bound for the A temperatures. 
The systematic error, caused by this variation was esti- 
mated to be 5%. The calculated feed-down contribution 
to the TT^ and p, p yields in Pb-|-Pb reactions is depicted 
in Fig. [6l 

The tracking inefficiency was estimated by simulation 
and reconstruction of embedded particles in real events, 
i.e. by determining the conditional probability of losing 
a track of given particle type and momentum, taking 
the proper track density environment into account. The 
inefficiency estimates are plotted in Fig. [T] These also 
include the inefRciency caused by the decay loss for the 
weakly decaying particles. 

The aim of our analysis is to measure the particle 
yields as a function of at y = 0, integrated over 
— 180° < 4> < 180°. Therefore, the particle spectra have 
to be corrected for the geometric acceptance. This cor- 
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FIG. 6: (Color online) Ratio of the estimated feed-down con- 
tribution to the raw measured particle yields. The top, center 
and bottom row of plots correspond, respectively, to Pb+Pb 
collisions with centrality ranges of 0-5%, 12.5-23.5% and 33.5- 
80% of the inelastic cross section (see Table |ll| . 



FIG. 7: (Color online) Detection inefficiencies for tt* 
(squares), (circles), p and p (triangles). The top, cen- 
ter and bottom row of plots correspond to Pb-|-Pb collisions 
with centrality ranges of 0-5%, 12.5-23.5% and 33.5-80% of 
the inelastic cross section (see Table |ll| . The decay loss for 
TT* and is included in the inefficiencies. 



rection can be performed without simulation, as the ac- 
cepted momentum space region is defined explicitly by 
our momentum space cut; only the (?/, (f>) shape of the 
particle distributions have to be known in each bin. 
The (f> distribution is uniform in each {y,Prp) slice due 
to the axial symmetry of particle production. For the 
y-dependence we use the fact, that the y spectra around 
midrapidity are approximately flat in Pb-|-Pb reactions 
(see [24]). Therefore, we assume a flat y distribution 
in —0.3 < y < 0.7 for the acceptance correction. The 
systematic error, caused by this approximation has been 
estimated to be below 2%. 

When all the discussed corrections are taken into ac- 
count, the Pb-fPb particle spectra bear the systematic 
errors listed in Table [TVl 



Points which correspond to a fitted amplitude of less 
than 25 particles or with a statistical error above 30% 
are not shown. Points below = 0.3GeV/c are not 
reported at all in order to avoid the ^ cross-over region. 



particle 


dE 
dx 


acceptance 


feed-down 


feed- down 


quadratic 


type 




correction 


yields 


shapes 


sum 


7r+ 


1% 


2% 






2.2% 


7r~ 


1% 


2% 






2.2% 


P 


1% 


2% 


3% 




3.7% 


P 


2% 


2% 


3% 


5% 


6.5% 


K+ 


4% 


2% 






4.5% 


R- 


1% 


2% 






2.2% 



TABLE IV: The list of estimated systematic errors and their 
sources on the values. 

dp 



IV. RESULTS 

A. Charged Hadron Spectra 

The resulting hadron spectra are shown in Fig. [51 
Preliminary versions, without the discussed corrections 
have been shown in [2^, [2^, [23| ■ The new spectra show 
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very good agreement with earlier low Pb+Pb results 
[2^ [28| . Also shown in Fig. [5] (lowest panel) are our 
7r+ , 7r~ spectra measured in p+p interactions [3] . 
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FIG. 8; (Color online) Invariant yields per inelastic collision 
of TT* (squares), K"^ (circles) andp,p (triangles) versus trans- 
verse momentum at midrapidity and ^/sj^ = 17.3 GeV 
energy. Results are shown for Pb+Pb collisions with central- 
ity ranges of 0-5% (top row), 12.5-23.5% (second row) and 
33.5-80% (third row) of the inelastic cross section (see Table 
lll|) . Previously published results for central collisions in the 
lower range [23 . [2^ are indicated by the bands in the top 
row. The bottom row shows published NA49 results [7t| for tv^ 
(left) and n~ (right) production in inelastic p+p collisions. 



B. Nuclear Modification Factors 

The nuclear modification effects are measured by the 
nuclear modification factors, which are particle yield ra- 
tios of two reactions, with appropriate scaling factors. 
The nuclear modification factors of a reaction Ai + Ai 
relative to A3 + A-^^ are most generally defined as: 



R 



Ai+A2/A3+Ai 



and 



^Ai+A2/Ai+Ai 



{N^^{A-i + Aj)) yield(v4i + A2) 
{N,^{A^+A2)) ■ yield(A3 + A4) ^ 



{Nw{Az + Ai)) yield(Ai + A2) 



{Nw{Ai + A2)) yield(A3 + A4) ' 



The ratio Ra+a/p+p is in the following abbreviated by 
the commonly used Raa^ while for R^^^/p+p the nota- 
tion RpA will be used. The scaling factors (N^^) and 
(Nw) are the average numbers of binary collisions and 
of wounded nucleons, respectively, which are calculated 
according to the Glauber model by Monte Carlo methods 
which simulate the geometrical configurations. The two 
different kinds of scaling assumptions are motivated by 
extreme scenarios of particle production. A pQCD-like 
particle production scheme would imply a binary collision 
scaling (particles are produced in binary parton-parton 
collisions). However, a scenario of a superposition of soft 
collisions (e.g. particle production via projectile excita- 
tion and decay) would rather suggest the scaling of par- 
ticle production with the number of wounded nucleons, 
as proposed in e.g. [2^ (each collision increases the en- 
ergy content of the excited nucleon) . As the particle pro- 
duction scheme has not been completely determined at 
SPS energies, especially below 2GeV/c transverse mo- 
mentum, both extreme scaling assumptions are exam- 
ined. 

To verify the presence or absence of nuclear effects, we 
calculate the modification factors Raa and Rcp (A+A 
central/peripheral) for hadrons using our particle spectra 
at y/s^ = 17.3 GeV. We compare our results to those 
obtained at y/s^ = 200 GeV, measured by the PHENIX 
experiment at RHIC 0, Q ■ As the production of 7r+ and 
7r~ is already very similar at our energy, we only show 
the average of tt^ for charged pions. The heavier parti- 
cle yields {p,p, , K^) proved not yet to be symmetric 
with respect to charge conjugation at our energy and we 
show them separately. The antiparticle/particle asym- 
metry of our spectra can be seen in Fig. [51 As p+p 
is a maximally isospin-asymmetric system, the 7r^/7r+ 
ratio is slightly below one. Pb-|-Pb is less asymmetric 
in isospin, thus its 7r~/7r+ production ratio is closer to 
unity. The isospin argument does not hold for the case 
of K~ / K'^ and p/p ratios, which are strongly infiuenced 
by the net-baryon density. 

The nuclear modification factor Raa is shown in Fig. 
[TOl for the tt* channel at top SPS and RHIC collision 
energies. The RHIC data are compared to RdA at the 
same center-of-mass energy, while for the SPS data RpAj 
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FIG. 9: (Color online) Ratios between antiparticle and par- 
ticle yields versus transverse momentum at midrapidity 
and ys„„ = 17.3 GeV energy. Results are shown for tv~ /tt^ 
(squares), K~ / K'^ (circles) and p/p (scaled up by a factor 
6, triangles). Plots correspond to Pb-|-Pb collisions with cen- 
trality ranges of 0-5% (top left), 12.5-23.5% (top right) and 
33.5-80% (bottom left) of the inelastic cross section (see Table 
ln| . Ratios from published NA49 results on 7r~ and tt"*" yields 
in p-(-p collisions 0] are plotted in the bottom right panel. 



measured at y/s^ = 19.4 GeV f9'|, is shown in addition. 
Since the RpA factors exhibit only a weak energy de- 
pendence for < < 4.5 GeV/c in the energy range 
19.4 < y/s^ < 27.4 GeV '^], the comparison to the 
.ys„„ = 17.3 GeV data is to a certain extent justified. 
However, the shrinking of the available momentum space 
region with decreasing collision energy makes this extrap- 
olation in energy somewhat uncertain, as the modifica- 
tion factors are expected to behave singularly close to the 
momentum space kinematic limit. 

The main observations with the binary collision scaling 
assumption are the following. Both the Raa and the 
RpA data show a similar low increase at y/s^ — 
17.3 GeV and = 200 GeV, up to p^ < 1 GeV/c. 

The RpA ratios rise above 1 at both collision energies with 
similar slopes, showing an excess of particle yield, which 
is often called Cronin effect. The 17.3 GeV Raa data 
keeps rising approximately linearly up to 2 GeV/c, where 
the statistics of the p+p reference spectrum runs out. In 
the case of 200 GeV Raa data, the ratio starts to decrease 
in this region, showing a large suppression both relative 
to the constant 1 line and to the RpA ratio (the "Cronin 
baseline"). The evolution of the 17.3 GeV data at > 
2GeV/c is not clear, because the available experiments 



FIG. 10: (Color online) Ratios of charged pion yields at 
midrapidity in central nucleus-|-nucleus and inelastic nu- 
cleon-|-nucleon collisions scaled by the number of binary col- 
lisions {R^^ , top row) and number of wounded nucleons 
{R^ , bottom row). The left column shows results from 
NA49 at the SPS {^/s;^ = 17.3 GeV) for the ratio of cen- 
tral (0-5%) Pb-|-Pb collisions to inelastic p-l-p reactions (filled 
squares) compared to the ratio of p-bW to p+p collisions at 
y/s^j^ = 19.4 GeV (the dotted line is shown to guide the 
eye). The right column shows measurements from RHIC 
i^A^ = 200 GeV) [2, H for the ratio of central (0-5%) 
Au+Au collisions (open squares) and d-|-Au reactions (open 
triangles). (The error bars attached to the constant 1 line 
indicate normalization uncertainty.) 



on p-l-p have no statistics beyond p^ — 2GeV/c. In the 
covered p^ region, the Raa ratio stays below the RpA 
Cronin baseline also at 17.3 GeV. 

When assuming wounded nucleon scaling, there are 
some differences compared to the above picture. Namely, 
the modification factors start from one at both energies, 
and show a rapid increase with similar slopes at the be- 
ginning. The RpA ratios stay below the Raa result as 
opposed to the binary collision case. The 200 GeV Raa 
data show a tendency of returning to the constant 1 line 
at « 3 GeV/c. 

A widely believed interpretation of the Cronin effect 
is momentum transfer from the longitudinal degrees of 
freedom to the transverse degrees of freedom by multi- 
ple scattering of partons or hadrons, depending on the 
picture. However, then this effect must also be present 
in A+A collisions, and may mask other nuclear effects. 
On the other hand the nuclear modification factor Rcp 
would not contain most of the contribution of the Cronin 
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effect. The nuclear modification factor Rcp is shown in 
Fig. [TT] for different particle types at top SPS and RHIC 
collision energies. When assuming binary collision scal- 
ing, the TT^ ratios show an amazing similarity at the two 
energies at < 2 GeV/c, above which the 17.3 GeV ra- 
tio seems to stay constant (below unity) , and the 200 GeV 
ratio begins to show a large suppression. The p ratios 
seem to saturate at > 2GeV/c at both energies and 
are almost identical over the whole range. The p ra- 
tios differ at the two energies. The K~^,K~ results are 
very similar at the two energies in the common p^ range. 
In fact, the NA49 results cover a slightly larger p^ 
interval than the published PHENIX measurements. 

When using wounded nucleon scaling, the modification 
factors start from one, and show a rapid increase. The 
TT^ , modification factors show differences between the 
top SPS and RHIC collision energies. The difference of p 
results is even larger. The p modification factors at top 
SPS and RHIC energies, however, are quite similar. 

In Fig. [T^ the available Rcp measurements for tt^ 
and 7r° at SPS and RHIC are compared [1,[E|33]. It is 
seen that the agreement between the result for tt^ from 
NA49 and 7r° from WA98 for p^ > 0.8 GeV/c is quite 
good within errors. The difference between the PHENIX 
TT^ and 7r° measurements is larger and outside errors. 

In conclusion, Rcp ratios for identified charged parti- 
cles scaled by the number of binary collisions are remark- 
ably similar at SPS and RHIC energies. Our present 
results at the SPS only reach the medium p^ region, 
where soft production processes and the quark coales- 
cence mechanism [sij may still dominate. It remains to 
be tested by future experiments, having access to an ex- 
tended p^ range, whether energy loss of hard scattered 
partons in dense matter is also a significant process at 
SPS energies. 



V. SUMMARY 

Invariant yields of tt^, p,p, particles were mea- 
sured in Pb-I-Pb reactions at ^/s^ = 17.3 GeV colhsion 
energy around midrapidity as a function of transverse 
momentum up to 4.5GeV/c, with overall systematic er- 
rors below fa 5%. Using these spectra and the previously 
published tt^ yields in p-l-p collisions, the nuclear mod- 
ification factor Raa for and Rcp for n^, p, p, 
were studied. 

The extracted Raa ratio, scaled with the number of 
binary collisions, shows a rapid rise with transverse mo- 
mentum in the covered region and does not show the 
strong suppression observed at RHIC. Interestingly, the 
Rcp ratios for 7r*,i^^, and p in the available region 
stay rather close to the y/s^ = 200 GeV RHIC results. 
In fact, Rcp ratios for p are almost identical. However, 
although the tt^ ratio stays below one also at the SPS, it 
shows a much larger suppression at RHIC above 2 GeV/c 
transverse momentum. When scaled with the number 
of wounded nucleons, the modification factors Raa and 



1 \ \ \ \ r~ 

Pb+Pb(0-5%) 
Pb+Pb{33.5-80%) 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 3.5 
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FIG. 11: (Color online) Ratio of yields in central to peripheral 
Pb-|-Pb and Au+Au collisions at ^s„„ = 17.3 GeV (squares) 
and 200 GeV 0] (triangles) . Left column shows ratio Rq'p 
using binary collision scaling, right column shows ratio B^p 
using wounded nucleon scaling. Centrality intervals are spec- 
ified in the legend of the top left panel. (The error bars at- 
tached to the constant 1 line indicate normalization uncer- 
tainty.) 
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Rcp start from one and show a rapid rise with p^, fur- 
thermore the RpA ratio stays below Raa in the covered 
region for both SPS and RHIC coUision energies. 



The range of our modification factor measurements 
is hmited by the statistics of the p+p and peripheral 
Pb+Pb collision reference spectra. A recently started 
second generation experiment, CERN-NA61 [32. IssL Isl] . 
will make it possible to extract reference spectra up to 
higher of 4 GeV/c in p+p and p+Pb collisions match- 
ing the range of the currently available central Pb-|-Pb 
data. These future measurements will allow a more sen- 
sitive test for the presence of high particle suppression 
at SPS energy. 
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